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ABSTRACT 

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radia- 
tive equilibrium with fixed surface fields between B — I0 12 and 2 x 10 13 G and effective temperatures 
logTcff = 5.5-6.8, as well as with surface B and T e # distributions around these values. The models 
arc based on the latest equation of state and opacity results for magnetized, partially ionized hydro- 
gen plasmas. The atmospheres directly determine the characteristics of thermal emission from the 
surface of neutron stars. We also incorporate these model spectra into XSPEC, under the model name 
NSMAX, thus allowing them to be used by the community to fit X-ray observations of neutron stars. 
Subject headings: radiative transfer - stars: atmospheres - stars: magnetic fields - stars: neutron - 
X-rays: stars 



1. INTRODUCTION 

Thermal radiation has been detected from ra- 
dio pulsars and radio-quiet neutro n stars (NSs ; 
see iKaspi. Roberts, fe Hardind l2006t lHaberll 120071 : 
Ivan Kerkwijk fc Kaplanll2007t IZavlinll2007l for reviews) 
and from soft gamma-ray repeaters and anomalous X- 
ray pulsars, which form the magnetar class of NSs en- 
dow ed with superstrong (B > 10 14 G) magnetic fields 
(see lWoods fc Thompsonll2006l for a review). Radiation 
from the surface of these NSs can provide invaluable in- 
formation on the physical properties and evolution of the 
NSs. Characteristics of the NS, such as the gravitational 
mass M, circumferential radius R, and surface temper- 
ature T, depend on the poorly constrained physics of 
the stellar interior, such as the nuclear equation of state 
(EOS) and quark and superfluid/superconducting prop- 
erties at supra- nuclear densities. Many NSs are also 
known to possess strong magnetic fields (B ~ 10 12 - 
10 13 G), with some well above the quantum critical value 
(B > Bq — 4.414 x 10 13 G). 

The observed radiation from a NS originates in a 
thin atmospheric layer (with scale height ~ 1 cm) 
that covers the stellar surface. To properly interpret 
the observations of NS surface emission and to pro- 
vide accurate constraints on the physical properties of 
NSs, it is important to understand in detail the ra- 
diative behavior of NS atmospheres in the presence of 
strong magnetic fields. The properties of the atmo- 
sphere, such as the chemical composition, EOS, and 
radiative opacities, directly determine the characteris- 
tics of the observed spectrum. While the surface com- 

1 Also at the Isaac Newton Institute of Chile (St. Petersburg 
Branch), St. Petersburg, Russia 



position of the NS is unknown, a great simplification 
arises due to the efficient g ravitational separation of 
light and heavy elements (seelAlcock fe Illarionovlll980l : 
iBrown. Bildsten. fe Chanel I2002D . A pure hydrogen at- 
mosphere is expected even if a small amount of fall- 
back/accretion occurs after NS formation; the total mass 
of hydrogen needed to form an optically thick atmosphere 
can be less than ~ 10 16 g. Alternatively, a helium atmo- 
sphere may be possible as a r esult of diffusive hydro- 
gen burning on the NS surface rtChang fe Bildsten! 120031 : 
IChang. Arras, fe Bildsten! [200l ) . Finally, a heavy ele- 
ment atmosphere may exist if no accretion takes place or 
if all the accreted matter is consumed by thermonuclear 
reactions. 

Steady prog ress has been made in modeling NS at- 
mospheres (s ee iPavlov et al]|1995l iHo fe Lai|[200li 120031: 
iZavlinl I2007L for more detailed discussion and refer- 
ences on NS atmosphere modeling). Since the NS 
surface emission is thermal in nature, it has been 
modeled at the lowest approximation with a black- 
body spectrum. Early works on atmospheric spec- 
tra assume emission from light element, unmagne- 
tized atmospheres (the latter assumption being valid 
for B < 10 9 G); computed spectra exhibit signifi- 
cant deviation from a Planckian shape and distinctive 
hardening with respect to a blackbodv dRomanil Il987t 
jRaiago pal & Romani 1996: Zavlin. Pavl ov, fc Shibanovl 
119961: Idansicke. Braje. fc Romani 2002f) . 

The strong magnetic fields present in NS atmospheres 
significantly increase the binding energi es of ato ms, 
molecules, and other bound states (see, e.g- JLafeOOlL for 
a review) . Abundances of these bound states can be ap- 
preciable in the atmospheres of col d NSs (i.e., those with 
surface temperature T < 10 6 K; lLai fc Salpeterl 119971 : 
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iPotekhin. Chabrier. fc S hibanov 199^). In addition, the 
presence of a magnetic field causes emission t o be 
anisotropi c (iPavlov et alj 1994t [Zavlin et al]|1995al) an d 
polarized (|Meszaros et al.l 119881: [Pavlov fc Zavlinll2000l) ; 
this must be taken into account self-consistently when de- 
veloping radiative transfer codes. The most comprehen- 
sive early studies of magnetic NS atmospheres focused 
on a fully ionized hydrogen plasma and moderate field 
stren gths (B ~ 10 12 -10 1 3 G;lMillerill99llShibanoy et all 
Il992t IPavlov et al.lll994h IZane. Turolla. fc Treves! l2000f ). 
These models are expected to be valid only for relatively 
high temperatures (T > a few x 10 6 K), where hydro- 
gen is almost completely ionized. More recently, atmo- 
sphere models in the ultra-strong field (B > 10 14 G) 
and relevant tempera ture r egimes have been presented 
(IHo fc Laill200lL[2^0l[Q^^cM[Zane et al.ll200HlLlovd 
2003 ; van Adelsberg fe Lail 200t: ; see iBezchastnov et al ' 
Bulik fc Miller! 119971 
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19961 iBulik fe Miller! 119971 for early work), and all of 
these rely on the assumpti on of a fully ion ized hydrogen 
composition (see, however. IHo et aLll2003l) . Mag netized 
non-hydro gen atmospheres have been stu died by iMillerl 
(|1992f > and lRaiagopal. Romani. fc Milled (|1997j ). but be- 
cause of the complexity of th e atomic physics, the mod- 
els were necessarily crude (see iMori lHlo1[2007l for more 
details). Only recently have self-consistent atmosphere 
models (|Ho et al.ll200"at IPotekhin et alJl2004t IMori fc Hd 
[200l using the l atest EOS and opacities for par tially 
ionized hydrogen (IPotekhin fc Chabriei)l2003l I2004D and 
mid-Z elements (jMori fc Hailevl 120021 f2006h been con- 
structed. 

Here we present a systematic tabulation of our 
partially ionized hydrogen atmosphere models for 
B = 10 12 -2 x 10 13 G. We incorporate these tables 
into XSPEC 2 (jArnaudl fl996l) . under the model name 
NSMAX, for use by the astronomical community. Note 
that the NS atmosphere models previously imple- 
mented in X S PEC are either non- magnetic (NSA G RAV : 
IZavlin et all fl99& NSSPEC: iGansicke et all 
NSATMOS: iMcClintock. Naravan. fc Rvbickil 
iHeinke eT al. 2006) or magneti c but fully ionized hydro- 
gen (NSA: IPavlov et al.l Il995f ): the last at two fields: 
B = 10 12 G and 10 13 G. In § El we give details on 
the construction of the atmosphere models. In § we 
present our results, and we summarize and mention 
future work in § |4j 

2. CONSTRUCTION OF ATMOSPHERE MODEL 

Thermal radiation from the surface of a NS is medi- 
ated by the stellar atmosphere. The model for emission 
through the atmosphere is constructed using a grid in 
Thomson depth r, photon energy E, and photon propa- 
gation direction (6k,<f>k), where 6k is the angle between 
the photon wave-vector k and the surface normal n and 
4>k is the azimuthal angle between k and the magnetic 
field B (see Fig. [TJ. Grid intervals are equally spaced 
logarithmically for depth 10~ 5 < r < 10 3 and energy 
0-01 < E < 10 keV and spaced every 5° for 6 k and 10° 
for <f>k (extra grid points are included around 9k = Ob, 
where Ob is the angle between n and £?); 6 grid points 
are used per decade in r, and 50 grid points are used 
per decade in E. Under typical conditions in NS at- 
mospheres with B > 10 12 G, radiation propagates in 

2 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 





Fig. 1. — Coordinate axes and angles used to describe the 
atmosphere model. Two Cartesian coordinate frames, (x n y n z n ) 
and (xyz), are considered. The z n -axis is directed along the surface 
normal n. B is the magnetic field, and ©s and <f>k are its polar 
and azimuthal angles in the (x n j/ n 2 n )-frame. fc is the photon wave- 
vector at the surface and lies in the (x n z n )-plane, Sj. is the angle 
between fc and n, and fc' is the photon wave-vector at infinity and 
is directed along z (fc = fc' in the absence of gravitational light- 
bending). 8 and tp are the polar and azimuthal angles of the surface 
point at the stellar surface in the (xj/z)-frame (clearly, the angle 
between n and fc' equals 8). For a magnetic dipole, rh is the unit 
vector along the magnetic axis, which lies in the (xz)-planc, 8 m is 
the angle between fc' and rh, and 7 is the magnetic colatitude of 
the surface point. 



two polarization modes (see, e . g. . iMeszarosll 1 9921 ) ; there- 
fore, the radiative transfer equations for the two cou- 
pled modes are solved to determine the emission prop- 
erties of a magnetic atmosphere. The self-consistency of 
the atmosphere model is determined by requiring that 
the deviations (at each Thomson depth) from radia- 
tive equilibrium and con stant total flux are C 1% and 
< 0.5%, respectively (see IHo fc Lall200ll IHo et al J 12003 
IPotekhin et al.ll2004L and references therein for details on 
the construction of the atmosphere models). The atmo- 
sphere models mainly depend on three parameters: the 
effective temperature T ff and the magnetic field strength 
B and inclination (see Fig. [IJ. The atmosphere 
models also have a dependence, through hydrostatic bal- 
ance, on the surface gravity g = (1 + z g ) GMj R 2 » 
1.328xl0 14 (l+z s ) (M/Mq) (R/W km)" 2 cm s~ 2 , where 
z g is the gravitational redshift, 1 + z g = (1 — r g /i?) -1 / 2 , 
and r g = 2GM/c 2 « 2.95 (M/Mq) km is the gravita- 
tional radius. Thus the atmosphere model depends on 
the NS mass M and radius R; however, the resulting 
spectra do not vary signi ficantly for differen t values of g 
around 2 x 10 14 cm s~ 2 (jPavlov et alJll995t see Figs. El 
E and©. 

The spectra from models constructed as discussed 
above only describe emission from either a local patch 
of the stellar surface with the particular T c s and 
B or a star with a uniform temperature and radial 
magnetic field of uniform strength at the surface. By 
taking into account surface magnetic field and tem- 
perature distributions, we can construct more physical 
models of emission from NSs; however these spectra 
from the whole NS surface are necessarily model- 
dependent, as the B and T distributions are generally 
unknown (see, e.g.. IZavlin et al.lll995aHZane et al.ll2001l: 
Ho fc Lai 2004: Per ez- Azorin et a l. 2006; Zane fc TuroTial 
Hoi 120071: see also iGeppert. Kiiker. fc Pagd 



lAguilera. Pons, fc MirallesT 



20071 



IPotekhin. Chabrier. fc Yako vlcv 2007; Re isenegger et all 
120071 and references therein, for recent work on surface 
B and T distributions). Therefore, we provide results 
from two sets of models: one set with a single B and T e s, 
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TABLE 1 

Magnetic Field Distribution of Neutron Star Surface 
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fully ionized hydrogen (see NSA; iPavlov et al.l Il995f) . 
and a set which is constructed with B and T c g varying 
across the surface. The latter is built by dividing the 
surface into regions with different B and T e ff- Relatively 
simple distributions are adopted (see Tables [1] and [2]) . 
In particular, we assume a dipolar magnetic field, after 
accounting for the effect of General Relativity, such that 
the surface distribution of B is given by 

B = (Spoie/2) [(2 + /)(«• rh)n - frh] (1) 

(iGinzburg fc Ozern oi 1965; see also IPavlov fc Zavlinl 
2000), where B po i c is the field strength at the magnetic 
pole, / = f(M/R) > 1 accounts for spatial curvature 
(/ = 1 in planar geometry), and rh is the direction of the 
magnetic moment. The surface te mperature distr i bution 
is calculated using the results of iPotekhin et all (2003) 
for the "canonical" NS mass M — 1.4M© and the radius 
R = 12 k m, corresponding to moderately stiff N S EOSs 
(see, e.g., lHaensel. Potekhin. fc Yakovlevll200"7[) . In or- 
der to minimize model dependence, we assume the T e g- 
distribution of an iron heat-blanketing envelope. This 
assumption does not change our results since, for any 
chemical composition of the envelope, the dependence 
of T e g on the magnetic colatitude 7 is similar to that 
given in lGreenstein fc Hartkd (|1983f) (see IPotekhin et all 
12001 . 

Emission from any point along a circle at a fixed 7 is 
given by the atmosphere model for that 7. Using equa- 
tion ([T]) and the coordinate transformation 

x n = — cos ip cos 9 x — sin ip cos 9 y + sin 9 z 

y u = — sin ip x — cos ip y (2) 

z n = cos ip sin 9 x + sin ip sin 9 y + cos 9 z , 

where 9 and ip are the polar and azimuthal angles of 
n with respect to the line of sight (see Fig. [T|) and 
(x, y, z) and (x n , y n , z n ) are the unit coordinate vectors 
for frames (x, y, z) and (x n , y n , z n ), respectively, one finds 
that 

cos 7 = cos ip sin 9 sin 9 m + cos 9 cos 9 m , (3) 
sin</>fc = sin 93 sin (9 m / sin 7, (4) 



B = B 



pole 



COS' 



: 7+(//2) 2 sin 2 7 , 



tan9s = (//2) tan7, 

where 9 m is the angle between 
4>k = <p in the special case 7 
ton wave-vector at infinity k' 



(5) 
(6) 

rh and k' and 
= 0. The pho- 
diffcrs from k at 



the surf ace due to gravitational redshift and light- 
bending dPechenick. Ftaclas. fc Coher]|1983l: |Pagill995l : 
IPavlov fc Zavlinll200Q[h " The latter is taken into account 
by ma king use of the approximation from iBeloborodovl 
(|2002f) (see also lZavlin. Shibanov. fc Pavlovlll995ri) . 



TABLE 2 

Effective Temperatures T e g/(10 5 K) of Neutron Star 
Atmosphere Models 
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Note. — Models assume g ■ 
for these temperatures use blackbodies with T 
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The spectrum for the entire surface is then computed by 
summing over the emission from different regions, 



F E = A 
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tt/2 



sin 9 k d9 k I 



kJ-EK 



;B,e B ), (8) 



cos (9 = (1 - cos0 fe )/(l - r s /R). 



(7) 



where Ie is the specific intensity. Note that E is the 
unredshiftcd photon energy. The explicit redshifting of 
the photon energy and flux spectrum is not done at this 
stage, though relativistic effects are taken into account 
[see, e.g., eqs. ((5]) and ((7j)]; redshifting is done in the 
XSPEC fitting code (see §[3])- At this point in the model 
calculation, the flux normalization A [= A(M,R,d), 
where d is the distance to the source] is taken to be unity 
(see §[3] for discussion of its model-dependence). 

We calculate Ie for four 7 values (see Table [1} and 
perform the integration in equation ([8]) by interpolat- 
ing between the calculated values. For T e ff(7 = 90°) < 
2 x 10 5 K, we use blackbody spectra at T = T e g (see Ta- 
ble [2]); the spectral contributions at these temperatures 
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contribute little to the total X-ray spectra [since T^g ^> 
?cff (7 = 90°); replacing the blackbody spectrum with 
zero values yields no appreciable difference in the result- 
ing integrated spectrum], which is dominated by emis- 
sion from the hotter regions of the NS surface. Strong 
absorption features, such as the proton cyclotron line at 
Esp = heB/rripC, are broadened due to the variation of 
B with 7. In order to reproduce this broadening in our 
interpolation, we first remap our calculated Ie as a func- 
tionof£/£ Bp : I E {6 kj fa; B, Q B ) = I(E/E Bp , 9 k , fa; 7). 
We then interpolate / in 9 k , fa, and B(j) for every fixed 
E/Eb p and substitute the resulting Ie into equation 

3. RESULTS 

For the first set of models (with uniform B and T e g), 
we consider g = 1.6 and 2.4 x 10 14 cm s -2 ; for a NS with 
M = 1.4M©, this corresponds to R = 12 km and 10 km, 
respectively. The magnetic field is B — 10 12 , 1.26 x 10 12 , 
2 x 10 12 , 4 x 10 12 , 7 x 10 12 , 10 13 , or 2 x 10 13 G and 
Ob = 0°. The effective temperatures span the range 
logT e ff = 5.5-6.8 (5.6-6.8 for B = 2 x 10 13 G) with the 
temperature interval between each model AlogT e ff « 
0.1. The temperature and abundance profiles for the 
atmosphere models with g = 2.4 x 10 14 cm s -2 and 
B = 10 12 G and 10 13 G are shown in Figures [2] and H 
respectively. The atomic fraction is the number of hydro- 
gen atoms with non-destro yed energy levels divided b y 
the total number of protons ( Potekhin & Chabrier 2003) . 
Here we only account for the ground-state atoms, which 
substantially reduces the computational work. This ap- 
proximation is justified because the fraction of atoms in 
excited states is small in most of the considered tempera- 
ture profiles: it does not exceed a few percent even when 
the abundance of ground-state atoms reaches tens of per- 
cent. The dependence of the atomic fraction on temper- 
ature and magnetic field is clear: lower temperatures or 
higher magnetic fields increase the abundance of bound 
species. The dependence on density is more complex: an 
increase in the atomic fraction with growing density (re- 
combination according to the modified Saha equation) 
competes with the decrease due to pressure ionization, 
which ultimately turns to complete ionization at high p. 

Figures HtiTDl show the (unredshifted) spectra for B = 
10 12 , 1.26 x 10 12 , 2 x 10 12 , 4 x 10 12 , 7 x 10 12 , 10 13 , 
and 2 x 10 13 G. The most prominent spectral fea- 
tures are due to the proton cyclotron line at Eb p = 
0.063 (B/10 13 G) keV, the s = -> 1 transition at 
E = 0.051 keV for B = 10 12 G and 0.14 keV for 10 13 G, 
the s = -> 2 transition at 0.075 keV for 10 12 G and 
0.23 keV for 10 13 G, and the peak of the bound-free tran- 
sition at 0.16 keV for 10 12 G and 0.31 keV for 10 13 G. 
Here s is the quantum number that measures transverse 
atomic excitations and corresponds to the projection of 
the angular momentum onto the magnetic field lines, 
whereas the longitudinal and Landau quantum numbers 
equal zero for the bound s tates involved in these tran- 
sitions (see IPotekhin! 1 19941 for a detailed description of 
the quantum numbers of a moving hydrogen atom) . All 
the spectral features due to atomic transitions are sub- 
stantially broa dened because of the "mo tional Stark ef- 
fect" (see, e.g.. iPotekhin fc Pavlovl [19971 and references 
therein). This "magnetic broadening" becomes stronger 
with increasing T and is another reason, in addition to 




Fig. 2. — Temperature and abundance profiles at various T c ff of a 
partially ionized hydrogen atmosphere with B = fO 12 G, 0s = 0°, 
and g = 2.4 X 10 1 * cm s~ 2 . 




Fig. 3. — Temperature and abundance profiles at various T c ff of a 
partially ionized hydrogen atmosphere with B = fO 13 G, 0s =0°, 
and g = 2.4 X 10 14 cm s" 2 . 

the decrease in the neutral fraction, for the disappear- 
ance of the features from the spectra at higher T e ff- 

For the second set of models, we take g — 1.6 x 
10 14 cm s -2 . The range of magnetic fields and effec- 
tive temperatures (-B,0B,T e ff) of the models are given 
in Table fJl These values correspond to a magnetic dipolc 
model of a NS with R = 12 km and M = 1.4M , in 
agreement with the chosen g. One subset of models has 
B = 10 12 G at the magnetic equator, while the other 
has B = 10 13 G at the pole. T ^ s is the mean effective 
temperature for the whole NS, which cor responds to the 
total heat flux from the surface (see, e.g. . IPotekhin et aLI 
2003). For each Tf$, the spectra are calculated as de- 
scribed in §[2] and using equation (JSj> . The resulting spec- 
tra are shown in Figures ITTHT41 for the cases 6 m = and 
90°. For comparison, the dotted lines show the atmo- 
sphere spectra for a uniform surface temperature and 
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Fig. 4. — Spectra at various T c ff of a partially ionized hydrogen 
atmosphere with B = 10 12 G and Sb =0°. Dotted and solid lines 
are for g = 1.6 and 2.4 X 10 14 cm s — 2 , respectively. 
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Fig. 5. — Spectra at various T c ft of a partially ionized hydrogen 
atmosphere with B = 1.26 X 10 12 G, S =0°, and g = 1.6 X 
10 14 cm s" 2 . 

radial magnetic field (see Figs. [4]and[9|). We see that 
the field distribution over the stellar surface substan- 
tially smears the spectral features from atomic transi- 
tions. This smearing is especially noticeable in FiguresfT3l 
andl!4[ where the atomic features are stronger due to the 
higher atomic fractions. 

Finally, we supply XSPEC, under the model name NS- 
MAX, with tables of the spectra shown in Figures H]- 
Q3J as well as a code to interpolate within each table; 
note that the model spectra with single (B, T G fr)-values 
span the photon energy range 0.05 < E < 10 keV, 
while the model spectra with (B, T g )-distributions cover 
0.09 < E < 5 keV. The code first unredshifts the energy 
bins of the observed spectrum, then obtains the fit spec- 
trum by linear interpolation via a weighted average of 
the nearest two model logT e ff and E, and finally red- 
shifts the fit spectrum by (1 + z g ) _1 . The code requires 




E (keV) 

Fig. 6. — Spectra at various T e ff of a partially ionized hy- 
drogen atmosphere with B = 2 X 10 12 G, Ob = 0°, and g = 
2.4 X 10 14 cm s" 2 . 




E (keV) 

Fig. 7. — Spectra at various T e g of a partially ionized hy- 
drogen atmosphere with B = 4 X 10 12 G, &b = 0°, and g = 
2.4 X 10 14 cm s" 2 . 

one switch parameter and two fit parameters (log T ff and 

1 + z g ); XSPEC automatically adds a third fit parame- 
ter (normalization A). The switch parameter indicates 
which table of model spectra to use: the differences be- 
ing due to the composition (only hydrogen at the present 
time; see §H]), B, Ob, and g for the first set of models 
and the composition, B, 6 m , and g for the second set 
of models. The normalization parameter A is conven- 
tionally taken to be equal to (R°°/d) 2 /(l + z g ) 3 , where 
R°° = R (1 + z g ), and the same R is used to calculate z g ; 
note that this prescription implies the emission region is 
the entire visible surface of the NS. 

4. SUMMARY 

We have constructed tables of model atmosphere spec- 
tra for neutron stars (with magnetic fields B = 10 12 - 

2 x 10 13 G and effective temperatures log Toff = 5.5- 
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Fig. 8. — Spectra at various T c g of a partially ionized hydrogen 
atmosphere with B = 7 X 10 12 G and S = 0°. Dotted and solid 
lines are for g = 1.6 and 2.4 X 10 14 cm s -2 , respectively. 
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Fig. 9. — Spectra at various T c g of a partially ionized hydrogen 
atmosphere with B = 10 13 G and 0g = 0° . Dotted and solid lines 
are for g = 1.6 and 2.4 X 10 14 cm s~ 2 , respectively. 

6.8) and incorporated these tables into XSPEC (un- 
der the model name NSMAX 3 ). These spectra are ob- 
tained using the most up-to-date equation of state and 
opacities for a partially ionized hydrogen plasma, and 
therefore they can describe emission from neutron stars 
with surface temperatures T < 10 6 K, where the abun- 
dance of bound species is appreciable, as well as neu- 
tron stars with T > 10 6 K. Thus we go beyond the 
previous magnetic neutron star spectral models provided 
in XSPEC, which assume fully ionized hydrogen atmo- 
spheres. Our implementation in XSPEC allows easy up- 
dates to the database of model spectra, so that tables of 
models with other magnetic field strengths (e.g., other 
than the seven fields, B = 10 12 , 1.26 x 10 12 , 2 x 10 12 , 
4 x 10 12 , 7 x 10 12 , 10 13 , and 2 x 10 13 G, currently pro- 




s' (keV) 

Fig. 10. — Spectra at various T e ff of a partially ionized 
hydrogen atmosphere with B = 2 X 10 13 G, 0s = 0°, and 
g = 2.4 X 10 14 cm s -2 . 




E (keV) 

Fig. 11. — Spectra at various total NS of a partially ion- 

ized hydrogen atmosphere covering a NS with R = 12 km and 
M = 1.4Mq, surface distributions of B and T c g according to the 
magnetic dipole model (see Tables [Tl and [2l Bp Q i e = 1.82 x 10 12 G 
and B = 10 12 G at the magnetic equator), and m = 0. Dotted 
lines correspond to atmosphere spectra with a uniform tempera- 
ture (T c ff = T^) and uniform radial magnetic field B = 10 G. 

vid ed) or other elem ents (e.g., carbon, oxygen, and neon; 
see lMori fc Holl2007l ) will be added as they become avail- 
able; the XSPEC-user merely specifies a switch parame- 
ter to indicate which set of models is to be used in the 
fitting. We have also constructed tables of model spectra 
that account for relativistic effects and dipolar magnetic 
field and temperature variations on the surface of the 
neutron star. These spectra are more realistic but also 
more model-dependent. They show significant smearing 
of spectral features compared to the models that assume 
a uniform magnetic field. 



3 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/nsmax.html We thank Keith Arnaud for assistance in incorporat- 
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ing NSMAX into XSPEC and the anonymous referee for 




helping to improve the clarity of the paper. WCGH 
appreciates the use of the computer facilities at the 
Kavli Institute for Particle Astrophysics and Cosmology. 
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nauka) grant NSh-2600.2008.2 and by RFBR grants 05- 
02-22003 and 08-02-00837. The work of GC is supported 
in part by CNRS grant PICS 3202. 
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Fig. 12.— The same as in Fig. [TTJ but with m = 90°. 
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Fig. 13.— The same as in Fig. [TT] but with S pole = 10 13 G 
and B = 5.5 X 10 12 G at the equator. Dotted lines correspond 
to atmosphere spectra with a uniform radial magnetic field B = 
10 13 G. 
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